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Section 1 


SUMMARY 

The work described in this report was carried out at Dynatech 

R/D Company for the National Aeronautics and Space Administration Lewis 

/ 

Research Center under Contract NAS-3-19716 with Mr. J. P. Joyce as the tech- 
nical monitor. 

The major purpose of the investigation was to select and measure 
the thermophysical properties of salts having potential for thermal energy 
storage to provide peaking energy in conventional electric utility power 
plants. 


Four candidate classes of electric power plants were considered. 

The systems and the selected salt for each were: 

1. Pressurized Water Reactor J (- liNO^ m.p. ~ 527K 

2. Boiling Water Reactor | | - 63 LiOH/37 LiCl eutectic m.p. ~ 533K 

3. Supercritical Steam Reactor - LiOH m.p. ~ 743K 

4. High Temperature Gas Reactor - Na2B^0^ m.p. ~ 1015K 

The thermal conductivity, specific heat (including latent heat 
of fusion), and density of each salt were measured for a temperature range 
of at least + 100K of the measured melt point. Measurements were made with 
both reagent and commercial grades of each salt. Where necessary techniques 
were developed to fabricate and prepare uncontaminated, close-to- theoretical- 
density solid samples for the measurements. The results are presented and 
discussed, in relation both to the very limited data available for any of 
the materials and to their potential use as thermal energy storage salts. 

Little difference in any property was found between the reagent and commercial 
grades of any of the four materials. 

The lithium salts are good candidate materials since they have 
relatively high specific heats and high latent heats of fusion (of the order 
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of 5 - 10 x 10 J/kg) . However, heat transfer Is poor by reason of the 
typically low values of thermal conductivity (of the order of 1 W/mK) . 

The thermal conductivity of the solid decreases with increase in temperature 
and there is an approximate 50% decrease at the melting point, followed by a 
small increase with increase of temperature. The density of the two lower 
temperature materials decreases at the melting point by over 20%. The lithium 
hydroxide decreases by only 2 %. 

Sodium diborate is a poor energy storage material since it does 
not appear to have a very high latent heat and has an apparent increase in 
density on melting. It behaves in an anomalous manner in the 500 to 700C 
regime, due presumably to the formation of a supercooled amorphous form. 

Following the measurements, a storage system was sized for association 
with a boiling water reactor of the class of the Braun SAR. The design conditions 
for the system were: 

Design Power 1250 MW 

Peaking Power 250 MW 

Storage Capacity (Electric) 600 MW-hr 

From the results of the measurements LiNO^ was preferred to the 
eutectic for this application. 

The design of this storage system required the following: 

Total mass of heat of storage material 19.2 x 10^ kg 

3 

Total volume of tank to contain heat 11600 m 

Total length of 0.025 m stainless steel tubing 1.7 x 10^ m 
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Section 2 


SELECTION OF HEAT STORAGE MATERIALS 

2 . 1 Concept of Storage System 

2.1.1 Power System Load 

Typically, electric power networks have a variable annual load 
demand, as shown in Figure 2.1 (Ref. 1). When the nature of this demand is 
such that excess base load generating capacity is available, energy storage 
systems, which can level some of the peaks and valleys of the time-variable 
load, become attractive. By more effectively using the most efficient generating 
equipment, energy storage systems can meet overall network electrical require- 
ments while reducing the cost of producing electricity and, possibly, the total 
energy input. 

Generating stations in an electric power network are classified 
as base load, intermediate cycling, or peaking, in accordance with the number 
of hours per year they are in operation. Base load plants are designed to be 
operated continuously throughout the year. In order to minimize their operating 
cost, base load plants are highly efficient, but relatively inflexible from 
the standpoint of power output. Intermediate cycling plants provide for daily 
cyclic loading, are generally somewhat less efficient and more costly to operate, 
and are brought on line between eight and twelve hours per day, or approximately 
3000 hours per year. Peaking plants are used only during periods of peak demand, 
are generally costly to operate, and are operated less than 1000 hours per year. 

2.1.2 Energy Storage 

Energy storage has been applied by electric utilities on a system 
basis for a number of years. Initially, when hydraulic power was an important 
component of the available system power, run of river storage was used to handle 
a large portion of the variable demand. This method was very effective, since 
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Typical Weekly Load Duration Curve 
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the initial cost of the facility was primarily in the dam and the efficiency 
of power production was relatively independent of rate. 

Pumped water storage is the principal method implemented in 
American electric utilities. Water is pumped into a reservoir during periods 
of low demand, or when available electric power is economically priced, and 
drawn out through a power turbine to provide peak loads. Pumped storage does 
introduce an inefficiency in the product of the component efficiencies of the 
pump and the turbine, normally between fifty and seventy percent. The total 
capacity available from this class of storage is limited by the natural sites 
available. 

A pumped storage unit is a separate component of the utility 
system. It receives and provides power to the electric network as a whole. 
Hence, its overall operating efficiency and cost should be penalized by the 
inefficiencies and costs of operating the electric network. 

2.1.3 High Temperature Thermal Storage 

High temperature thermal storage, the alternate technique under 
consideration in this study, receives energy in the form of heat rather than 
electricity. Therefore, it is intimately related to a specific power plant. 

Thermal storage receives excess heat generated in a relatively 
efficient unit and gives this energy to a power conversion system during 
periods of peak demand. Because of their high efficiency and low flexibility, 
base load plants with periodic excess capacity are ideal locations for high 
temperature thermal storage systems. To the extent that excess capacity 
exists in base load plants in a network, thermal energy storage systems can 
augment or replace less efficient, more costly peaking equipment. 

Thermal storage could be applied to a cycling power plant. On 
a twelve hour cycling plant, heat would go into storage half the time and be 
drawn from storage and the heat generator during the other half. The power 
producer in this case would be twice the size of the power producer for a 
base load power plant using the same heat source. 
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The design conditions for this case are specified by the project 
monitor for a peaking application (Reference 2). These recommendations are 
given in Table 2.1. 


Table 2.1 


STORAGE 

SYSTEM DESIGN 

SPECIFICATIONS 


Absolute 

% of Daily Base Load 

Base Load Plant, MW 

e 

1000 

100% 

Peaking Power Rate, MW g 

200 

20% 

Storage Capacity, MW^-hr 



(Electric) 

480 

2% 


The two percent of equivalent electric output to storage requires 
more than two percent of the heat input. If the storage is being filled continually 
by an increased production of thermal energy by the heat source, this increase 
amounts to about three percent of plant design capacity (two percent divided 
by the estimated efficiency of the overall peaking system relative to the 
base load plant). If the storage is being filled during five hours of reduced 
demand, the reduction in plant output is twenty to thirty percent. The heat 
input to the cycle is reduced by 15% as the storage unit is building up its 
ability to produce electric power at the rate of 10% of base load capacity. 

In addition, there is a reduction in the cycle efficiency at part load. 

2.1.4 Classes of Storage Materials 

High temperature thermal storage can use either sensible heat 
or heat of fusion as an energy receiver. Heat of fusion has the advantage 
of permitting a relatively small receiver design, since a large amount of 
energy can be stored per unit volume of the storage material. It has the 
disadvantage, as compared to a liquid storage medium, of requiring heat transfer 
through a solid which is likely to be a poor conductor. In this contract we are 
considering only heat of fusion materials. 
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Candidate materials include pure substances and eutectic mixtures 
which solidify at a constant temperature. This permits design at a maximum 
effective temperature difference between the source temperature and the tem- 
perature at which the peaking power system receives heat, as shown in Figure 2.2. 
Non-eutectic mixtures may have a narrow melting temperature which would not 
affect the temperature difference appreciably. However, during melting and 
solidification there is a gradual segregation of the mixture components, 
since the equilibrium mixture ratios are not equal at the solid-liquid inter- 
face. One component is preferentially depleted from the solid, and the char- 
acteristics change over time. Thus non-eutectic mixtures are not desirable. 

2 . 2 Factors Affecting Storage Material Selection 

High temperature thermal storage is a component of a system to 
provide peaking power for an electric utility network. The ultimate decision 
on the application of this concept is economic, based on the total cost of 
such a system as compared to a comparable peaking plant to satisfy the same 
region of the power spectrum, and is conservative, based on possible savings 
in oil and gas typically used in peaking plants. The final selection of the 
material is dependent upon its effect on the cost of the storage component. 

In the process of material selection, without a specific design under con- 
sideration, it is only possible to determine the trends of this effect, in 
terms of relative values for different candidate substances. 

The three components of the thermal storage system are: the heat 

transfer component that transfers heat from the prime source, the heat transfer 
component that transfers heat to the power cycle, and the storage vessel that 
contains the heat storage material. The precise relationship between the cost 
of these components depends upon the method of design integration. However, 
the cost of the two heat transfer units generally is proportional to heat 
transfer area, and the cost of the storage vessel is an increasing function 
of total volume. The cost of the storage material is directly proportional 
to the required volume. 
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The principal property of the storage material which determines 

the heat transfer area is the thermal conductivity of the solid. Heat is 

transferred from a flowing cycle fluid through the vessel, through a solid 

layer of the storage material, to the solid-liquid interface. The same process 

occurs in reverse when the heat is supplied to the peaking cycle. Therefore, 

heat has to be transferred through an accumulating or decreasing layer of 

solid. To the extent that this layer controls the heat transfer, the desirable 

1/2 

heat transfer parameter is (Ah^p) (see Appendix A). If the storage vessel is 

separate, the cost of the containment vessel is a function of a volume parameter, 

h p. If the heat transfer surface acts as the outside wall of the storage 
r 

vessel, the total cost of material required determines the storage cost para- 
meter. This is proportional to h /C. 

r 

Other properties also affect cost. Corrosiveness of the salt 
determines the compatible materials and thus the costs. The expansion during 
melting determines the necessary expansion provisions to be incorporated into 
the design. The specific heat may affect the useful enthalpy change during 
melting if a wide temperature range is used. 

With the limited selection of materials in the proper temperature 
ranges, and limited data availability on thermal conductivity, the cost para- 
meter and the volume parameter are used for selection. Excessive corrosiveness 
is cause for elimination, as is lack of availability in large quantities. 

2.3 Systems to be Considered 

Four classes of power plants are candidates for high temperature 
thermal storage systems. These are: 

1. Boiling water nuclear reactors. 

2. Pressurized water nuclear reactors. 

3. Supercritical, fossil fuel fired, steam systems. 

4. High temperature gas reactors. 


9 



This selection of candidates leads to three temperature ranges 
of interest. These are: 

Water-cooled reactors 910 - 980R 

Super critical steam 1410 - 1535R 

Gas cooled reactors 1760 - 1860R 

Their thermodynamic cycles determine the appropriate temperature 
range for the given application. 

The precise temperature range depends upon the cycle configuration. 
Figure 2.3 shows two potential overall configurations: one with the storage 

system in parallel; the other, in series. From the standpoint of the power 
plant, it is desirable to have the storage system in parallel because it 
minimizes the effect on the power plant system. From the standpoint of the 
storage system, it is better to have the series arrangement if there is an 
appreciable temperature drop in cycle working fluid. (A series arrangement 
with the storage unit following the power system puts a maximum penalty on 
the storage system.) Table 2.2 gives the minimum temperature of the fluid 
supplied to the storage system for each mode of operation described above. 

Since the water cooled reactors supply heat over a small temperature range, 
the mode of operation has little effect on the design temperature. For the 
higher temperature systems, parallel operation severely limits the temperature 
range of operation. 

The pressurized water reactor and the boiling water reactor have 
virtually the same temperature requirements. Two materials are selected, 
either of which may prove to be superior when more definitive information 
about properties is available. 
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Figure 2 . 3 


Cycle Configurations 
Pressurized Water Reactor 



Parallel 
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Table 2.2 


STORAGE DESIGN TEMPERATURE RANGE FOR SPECIFIC POWER PLANTS 



Boiling Water 
Reactor 

Pressurized 
Water Reactor 

Supercritical 

Steam 

High Temperature 
Gas Reactor 

Contractual Temperature Range 
°K 

High 

539 

544 

853 

1033 

Low 

506 

511 

783 

978 

Temperature Range 
of Working Fluid °k 

High 

559 

591 

811 

1040 

Low 

551 

563 

551 

603 

Minimum Temperature Available 
to Storage for Several 
Cycle Configurations 

Series, Constant 
Heat Removal 

559 

590 

803 

1027 

Series, Nighttime 
Heat Removal 

558 

587 

772 

974 

Parallel 

551 

563 

551 

603 


































2.4 Comparison of Potential Heat Storage Materials 


2.4.1 Materials Under Consideration 

A review is made of all potential heat of fusion materials over 
the range of temperatures from 450 K to 1150 K. Table 2.3 lists these materials, 
the corresponding melt temperatures, estimated heats of fusion, and densities. 
This list exceeds the potential operating ranges of the cycles under consider- 
ation and is not limited to the temperature range for the specific cycles. 

This is a reference list given for future application of heat storage materials 
at other temperatures. Materials omitted from this list are not considered. 

Table 2.4 gives additional properties, where available, for the 
materials in Table 2.3. Since there is only a limited amount of thermal con- 
ductivity data available, it is difficult to apply the heat transfer parameter 
for comparison. Thus only the volumetric heat of fusion, which is a factor 
in this parameter, is used. 

2.4.2 Selection of Materials for Water Cooled Reactors 

Table 2.5 lists the candidate materials in the specified temper- 
ature range and gives the effective latent heat on a cost and volume basis. 

Sodium hydroxide has too high a melt temperature, and the next three candidates 
are marginal in melt temperature, being incompatible with the boiling water 
reactor. All materials are available in car load lots, and none of the materials 
exhibit major compatibility problems. 

NaNO^-NaCl among the three marginal candidates is eliminated 
because its volumetric heat of fusion is low. The other two are comparable, 
both on a cost and volume basis. 

LiOH-NaOH is eliminated because it has a low volumetric latent 
heat, a low temperature of fusion, and a high relative cost. Both LiNO^ and 
LiCl-LiOH are expensive, but they have good volumetric heat of fusion properties. 
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Table 2.3 

HIGH TEMPERATURE HEAT OF FUSION MATERIALS 


No. 

Material 

Melt Temp. 
°K 


Density 

kg/m 3 

1 

KF 

1125 

454 

2480 

2 

Na 2 C °3 

1125 

279 

2530 

3 

Ca 

1123 

221 

1540 

4 


1121 

1044 

2640 

5 

LiB0 2 

1108 

698 

1400 

6 

75 NaF + 25 MgF 2 

1105 

649 

2680 

7 

62.5 NaF + 22.5 MgF 2 + 15 KF 

1082 

607 

2630 

8 

NaCl 

1074 

484 

2180 

9 

Cal 2 

1057 

142 

3490 

10 

CaCl 2 

1046 

256 

2280 

11 

kci 

1043 

372 

1990 

12 

67 LiF + 33 MgF 2 

1019 

947 

2630 

13 

65 NaF + 23 CaF 2 + 23 MgF 2 

1018 

574 

2760 

14 

Na 2 B 4°7 

1013 

523 

2370 | 

KS1 


998 

605 

2200 

16 

MgCl 2 

988 

454 

2240 

tsm 

60 KF + 40 NaF 

983 

479 

2510 

lie 

LiH 

956 

2582 

790 

urn 

A1 

933 

388 

2710 

Em 

60 LiF + 40 NaF 

925 

816 

2480 

tsm 

Mg 

923 

372 

1740 

22 

46 LiF + 44 NaF + 10 MgF 2 

905 

858 

- 

2610 

23 

52 LiF + 35 NaF + 13 CaF 2 

888 

640 

2630 

warn 

LiCl 

883 

470 

2070 

25 

52 NaCl + 48 NiCl 2 

843 

558 

2840 


26 

Ca(N0 3 ) 2 

834 

130 

2500 

mam 

73 LiCl + 27 NaCl 

825 

430 

2090 

28 

48 NaCl + 52 CaCl 2 

773 

328 

2160 

wsm 

49 KF + 51 LiF 

765 

461 

2560 

j 30 

1 i 

80 Li 2 C0 3 + 20 K 2 C0 3 

763 

377 

2170 

E9B 

LiOH 

743 

930 

1460 

Em 

11.5 NaF + 42 KF + 46.5 LiF 

727 

442 

2560 

33 

NaCl + MgCl 2 

713 

326 

2240 

EB 

80 LiOH + 20 LiF 

700 

1163 

1550 

35 

KOH 

673 

140 

2040 

36 

LiCl + KCI 

623 

255 

2030 

37 

kno 3 

613 

128 

2110 

Em 

NaOH 

593 

160 

2070 

39 

S3S99HBSMHHHHHNS! 

588 

244 

1730 

40 

93.6 NaN0 3 + 6.4 NaCl 

568 

191 

2260 

41 

95.3 NaOH +4.7 Na„S0, 
2 4 

566 

326 

- 

42 

7.8 NaCl + 6.4 Na 2 C0 3 + 85.8 NaOH 

555 

316 

2100 

m 

37 LiCl + 63 LiOH 

535 

437 

1640 

44 

NaCl + ZnCl 2 

533 

198 

2480 

warn 

23 LiOH + 77 NaOH 

528 

233 

1890 

46 

LiN0 3 

527 

379 

2400 

47 

A1C1 3 

468 

290 

2440 

E9 

NaOH + KOH 

463 

233 


49 

Li 

453 

442 

530 





























































































































































































Table 2.4 


AVAILABLE THERMAL PROPERTIES OF HEAT STORAGE MATERIALS 


No. 

Material 

Specific Heat 
kj/kgK 

Thermal Conductivity 
W/mK 

46 

LINO 3 

1.63 

- 

38 

NaOH 

1.47 

1.54 

26 

Ca(N 0 3 ) 2 

0.88 

- 

19 

A1 

0.92 

204.2 

18 

LiH 

8.04 

7.0 

16 

MgCl 2 

0.75 

- 

15 

L 1 2 C °3 

2.64 


13 

65 NaF + 23 CaF 2 + 12 MgF 2 

1.17 


12 

67 LiF + 33 MgF 2 

1.42 

- 

43 

37 LiCl + 63 LiOH 

0.75 


42 

7.8 NaCl + 6.4 Na 2 C0 3 + 85.8 NaOH 

2.51 

- 

34 

80 LiOH + 20 LiF 

2.14 

2.6 


15 






































Table 2.5 


MATERIALS CONSIDERED FOR WATER COOLED REACTORS 


No. 

Material 

Melt Temp. 
°K 

Cost Factor 
kJ/$ 

Volume Factor 
kJ/1 

m 

NaOH 

593 

992 

965 

2 

93.6 NaN0 3 +6.4 NaCl 

568 

1245 

435 

3 

95.3 NaOH +4.7 Na„SO, 

2 4 

566 

928 

703 

■ 

7.8 NaCl +6.4 Na 2 C0 3 +85.8 NaOH 

555 

1234 

674 

5 

37 LiCl + 63 LiOH 

535 

100 

728 

6 

NaCl + ZnCl 2 

533 

492 

469 

m 

23 LiOH +77 NaOH 

528 

177 

456 

8 

LiN0 3 

527 

111 

883 







































whereas the less expensive chloride eutectic has a low volumetric heat of 
fusion. Since the volumetric heat of fusion is considered the more important 
parameter, the selected pair is: 

1. lino 3 

2. 37 LiCl +63 LiOH 

However, the entire group of materials, with the exception of the two at the 
temperature limits, would be satisfactory candidates if considered only for 
the pressurized water reactor. Those with low volumetric heats of fusion have 
low costs which tend to offset the disadvantage. 

2.4.3 Selection of Materials for Supercritical Steam System 

Table 2.6 shows the considered candidate materials. The first 
three are present only because they fall within the originally specified tem- 
perature range, which proved to be too high. The last four were not evaluated 
in the original selection process. The chlorides have very low volumetric heats 
of fusion. The five f-luoriHes-*re competitive in performance with the lithium 
hydroxide, while sharing the two disadvantages of being eutectics and contain- 
ing highly corrosive components. 

In the original evaluation, the decision between the carbonate and 
the lithium compounds depends almost completely on their relative corrosion 
characteristics, since the heat of fusion of the carbonate is lower by a 
factor of two. From discussions with Mr. Gary Drage of the Naval Undersea 
Warfare Center, Mr. LeRoy Grantham of Atomics International and Mr. Worth 
Percival, Private Consultant, 707 Barham Down, Manchester, MO. 63011 (314- 
227-7923) , it was determined that the lithium compounds do not present in- 
surmountable containment problems at fusion temperatures. Reference 3 gives 
data on the eutectic, demonstrating that nickel corrosion rates are low enough 
to permit a twenty year life. Mr. Percival, citing his experience with LiOH, 
states that failures were caused by extraneous factors. Therefore the carbonate 
is eliminated. 
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Table 2.6 


MATERIALS CONSIDERED FOR SUPERCRITICAL STEAM SYSTEMS 


No. 

Material 

Melt Temp. 
K 

Cost Factor 
kJ/$ 

Volume Factor 
kJ/1 

1 

52 NaCl + 48 NiCl 2 

843 

33 

1585 

2 

Ca(N0 3 ) 2 

834 

- 

325 

3 

73 LiCl + 27 NaCl 

825 

206 

899 

4 

80 Li 2 CX> 3 + 20 K 2 C0 3 

763 

13 

818 

5 

LiOH 

743 

199 

1358 

6 

80 LiOH + 20 LiF 

700 

203 

1803 

7 

48 NaCl + 52CaCl 2 

773 

- 

708 

8 

49 KF + 51 LiF 

765 

- 

1180 

9 

11.5 NaF + 42 KF + 46.5 LiF 

727 

- 

1132 

10 

NaCl + MgCl 2 

713 

- 

730 







































The two advantages of the eutectic are the higher heat of fusion 
and the negligible volume change. The thermodynamic test results given in 
Reference 4 are questioned for three reasons: 

1. Mr. Drage has found the conductivity measurements made by 
the same contractor are wrong. 

2. The Dynatech Thermatest Department has not encountered a 
material for which the specific heat of the solid exceeds 
that of the liquid. 

3. In general, the heat of fusion of a eutectic is less than that 
of either component. 

Thus the decision is based on the relative importance of the 
design difficulties raised by the negative expansion versus the potential 
problems of handling a eutectic, where precise control of composition is 
critical. Since Mr. Percival has successfully designed and operated LiOH systems 
for a similar application over a long period of development, LiOH becomes the 
recommended substance. 

2.4.4 Selection of Material for Gas Cooled Reactor 

Table 2.7 shows the comparison among the candidate materials for 
use with the Gas Cooled Reactor. One of the group, LiCO^, can be eliminated 
because of relative expense. One, MgCl, can be eliminated because it has a 
low volume factor with no redeeming features. Satisfactory designs can be 
based around any of the others . 

Two, LiH and 67 LiF + 33 MgF^, have very high volumetric heats of 
fusion. The low volume factor for aluminum is offset by its high conductivity, 
which is at least one order of magnitude higher than any of the others. That 
conductivity is so high that it is doubtful that the system design will be 
limited by heat transfer through the solid layer. There is extensive previous 
work to determine the properties of aluminum (Refs . 7, 8) and of lithium hydride 
(Refs. 9, 10, 11). Repetition is not considered to be desirable. 
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Table 2.7 


MATERIALS CONSIDERED FOR HIGH TEMPERATURE GAS REACTORS 


No. 

Material 

— 
Melt Temp. 

K 

— 
Cost Factor 

kJ/$ 

Volume Factor 
kJ/1 

1 

67 LiF + 33 MgF 2 

1019 

237 

2489 

2 

65 NaF + 23 CaF 2 + 12 MgF 2 

1018 

883 

1582 

3 

Na 2 B 4°7 

1013 

2560 

1214 

4 

Li 2 C °3 

998 

17 

1328 

5 

MgCl 2 

988 

347 

1012 

6 

LiH 

956 

123 

2028 


A1 

933 

388 

1051 














The two fluoride eutectics and sodium borate give a range of 
potential properties. As the expense of the storage material increases, the 
volumetric heat of fusion goes down. The final decision is based on other 
considerations. Sodium borate is known to be an inert material and is readily 
available, as indicated by the low cost, whereas the fluorides are noted for 
difficulty in handling and corrosiveness. Therefore, the borate is selected. 


The principal difficulty expected with the borate arises from its 
tendency to form a glassy solid, sub cooling as heat is removed. This subcooling 


may be as much as 200°F unless a nucleating agent is developed for this material. 
Thus it is essential that a study of nucleation proceed early in the development 
test phase for the thermal storage system. 


2.4.5 Summary of Selection 


The materials selected for properties testing are: 


jlication 


Melt Temperature, 


Pressurized Water Reactor ) 

LiN0 3 

527 

Boiling Water Reactor ) 

37 LiCl + 63 LiOH 

533 

Supercritical Steam 

LiOH 

743 

High Temperature Gas Reactor 

Na 2 B 4°7 

1015 
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Section 3 


MEASUREMENT OF THERMOPHYSICAL PROPERTIES 

3.1 Outline of Investigation 

The necessary properties to be evaluated in the overall study 

were: 

1. melt temperature or range, 

2. thermal conductivity in solid and molten phases, 

3. heat capacity in solid and molten phases including latent 
heat of fusion, and 

4. density and coefficient of volumetric expansion in solid and 
molten phases. 

For each material both reagent and commercial grades were to be evaluated 
over an approximate range of 100K above and below the melting points of each. 

The overall temperature range for the four materials was between 500K and 
1050K depending upon the application of a particular material. 

Standard accepted techniques were chosen for the measurement of 
each property, with modifications where necessary due to associated sample- 
related problems. Of the properties desired, the thermal conductivity, 
especially in the molten phase, was the one which presented the major problems 
and where potential errors would be greatest. 

In general, the chosen salts were highly reactive both in air and 
with many materials. Thus special handling of materials and samples was 
necessary at all stages of production, assembly and measurement and in some 
cases special materials and techniques were required both in the fabrication 
of test samples and in the subsequent measurements. This latter point is not 
covered in detail in Section 3.3 dealing with the sample preparation procedures. 
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The individual methods or techniques chosen were : 


1. Measurement of melt temperature or range by differential 
thermal analysis, in order to determine both the overall tem- 
perature range of testing and how a particular measurement 
was to be carried out. 

2. Measurement of thermal conductivity in the solid phase by the 
comparative flat slab method and in the molten phase by a 
special modified comparative technique. 

3. Measurement of specific heat, including latent heat of fusion, 
by adiabatic calorimetry for the two materials being studied 
in the temperature range 400 to 600K and by drop calorimetry 
for the other two higher temperature materials. 

4. Measurements of density and coefficient of expansion in the 
solid state by a combination of dilatometry and pycnometry 
and in the molten phase by a buoyancy technique. 

Measurements of melting point, molten density and of specific heat could be 
determined on the materials in the crystal or powder form. However for the 
solid and molten thermal conductivity and solid linear expansion determinations 
measurements had to be carried out on solid crystalline samples prepared from 
the original crystals or powders. Details are given in Section 3.3. 

3.2 Details of Materials 


As a result of the initial analysis described earlier the fol- 
lowing salts were evaluated: 

1. lithium nitrate 

2. the 63% lithium hydroxide, 37% lithium chloride eutectic composition 

3. lithium hydroxide 

4. sodium tetraborate 


23 



The raw materials were obtained in sealed containers from various 
suppliers. A chemical analysis of each lot was also obtained prior to any 
measurements. Details of the sources and of the compositions of each are 
given in Tables 3.1 and 3.2 respectively. 


3.3 Preparation of Specimens 


met: 


In the preparations of the specimens, two conditions had to be 


1. Fabricate a solid, 100% dense specimen with the proper geometry. 

2. Prevent contamination of the sample material. 

A survey of available materials to contain the molten salts was made to find a 
suitable one for each. Besides the noble metals, where costs are prohibitive, 
only nickel, high-nickel alloys, and some stainless steels had acceptable cor- 
rosion rates. Since molten salts will also corrode metal oxides more readily 
than the base material, the candidate material had to be resistant to oxidation. 

Type 304 stainless steel was chosen finally for preparation of 
all materials. It satisfied the basic requirements and was readily available 
in the stock sizes required for preparation of the appropriate thermal con- 
ductivity and linear expansion samples. However during the course of the 
preparation of the sodium borate considerable difficulty was experienced with 
obtaining the release of the sample from the mold. Following various unsatis- 
factory experiences with the material, ATJ sulphur-free graphite had to be used 
for this one material. A special quartz glass had also been tried but unsuccessfully. 

The thermal conductivity samples for measurement in the solid 
state were approximately 51 mm diameter and 10 to 12 mm uniformly thick with 
a 0.3 mm groove cut across each smooth flat surface. For the molten state 
the sample was approximately 75 mm diameter and 20 mm thick. The linear ex- 
pansion samples were 6 mm diameter and 50 mm long and the ends were machined 
plane parallel. 
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Table 3.1 


DETAILS OF SOURCES OF MATERIALS EVALUATED 


Material 

Grade 

Supplier 

Lot No. 

Physical Description 

Lithium Nitrate 

Reagent 

J. T. Baker 
Chemical Company 

35815 

White Crystalline 

Lithium Nitrate 

Commercial 

Lithium Corporation 
of America 

419-101 

White Crystalline 
Bulk Density 1440 kg m 

Lithium Hydroxide 

Reagent 

Mathes on , Coleman 
& Bell 

B 8F25 

White Crystalline 

Lithium Hydroxide 

Commercial 

Lithium Corporation 
of America 

406-27-13 

White Anhydrous Powder 
Bulk Density 450 kg m 

Lithium Chloride 

Reagent 

J . T . Baker 
Chemical Company 

405688 

White Crystalline 

Lithium Chloride 

Commercial 

Lithium Corporation 
of America 

411-66 

White Anhydrous Crystalline 
Bulk Density 995 kg m 

Sodium Tetraborate 

Reagent 

. J. T. Baker 
Chemical Company 

508555 

White Crystalline 

Sodium Tetraborate 

Commercial 

U.S. Borax 
Company 

N/A 

White Crystalline 



K) 

ON 


Table 3.2 

CHEMICAL COMPOSITION OF THE "AS RECEIVED" MATERIAL 








































































Molds with the required dimensions were then fabricated. They 
were designed in three parts. A flat plate with a polished surface was used 
as a base. On this sat a cylinder of an appropriate inner diameter with the 
inner surface polished and tapered to simplify the removal of the sample slug. 

A clamping device sealed the tube and plate together. Once the specimen was 
cast, the base plate could be removed and the specimen pressed out of the 
containing cylinder. 

To prepare a sample the following procedure was finally chosen. 

The specimen material was dried to constant weight at HOC in vacuum. Correct 
proportions of material were weighed out and placed in a nickel beaker and 
then placed into a sealed oven purged with an inert gas. The oven was then 
heated slowly to above the appropriate melt temperature of the particular salt. 

As the specimen material was being melted, the mold was heated to approximately 
the same temperature as the salt. The molten salt was poured carefully into 
the mold. The base of the mold was cooled such that crystallization occurred 
only on the bottom surface. The top surface was held above the melting point 
by a heater placed on top of the mold to create a temperature gradient through 
the salt. By slowly reducing the temperature, the rate of crystallization could 
be controlled. As soon as a suitable thickness had been crystallized, the entire 
assembly was allowed to cool to room temperature. The mold was then disassembled 
and the specimen pressed out, ready for subsequent surface preparation. 

This procedure produced samples very close to theoretical 
maximum density with very few minute voids. All other techniques provided too 
many voids within a sample. Once a test material was taken out of its disposable 
sealed container it was handled inside plastic glove boxes in an inert atmosphere 
of dry argon during all phases of sample preparation, fabrication, and assembly. 
Additional details relative to each salt are given as follows. 

3.3.1 Lithium Nitrate 

The low melting temperature and corrosion rate made lithium nitrate 
the easiest material to mold. Large differences in density from the solid to 
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the molten state made it difficult to obtain a specimen free of voids. By 
carefully controlling the rate of crystallization, specimens 92 to 95% dense 
were finally fabricated. 

3.3.2 63 LiOH / 37 LiCl Eutectic 


This was similar to the lithium nitrate in that large differences 
in the density of the two phases produced voids in the cast slug. Controlling 
the rate of crystallization improved this situation. Each component was weighed 
in correct amounts and melted individually. They were poured together and 
mixed and allowed to solidify. The slug was crushed into a powder and then 
remelted in a new cup and then normal preparation procedure followed. The 
finished samples had a very slight greenish tinge, due presumably to oxidation 
reaction with the nickel in the case of melting the lithium chloride. The 
material was very brittle and during the course of subsequent final sample 
machinery and preparation a number of samples cracked before suitable test 
specimens were finally made available. 

3.3.3 Lithium Hydroxide 

Small differences in density between the two phases simplified 
the casting of dense specimens. However, corrosion of the mold made removal 
of the sample slug very difficult. The slug bound itself to the pitted sur- 
face of the mold and shattered when it was pressed out. The final solution 
to the problem was to machine and polish all contacting surfaces each time they 
were used and minimize the time that the salt remained molten. 

3.3.4 Sodium Tetraborate 

This presented the most serious problems in the fabrication of 
solid 100% dense crystalline specimens. It was found to subcool to form a 
glassy, non-crystalline mass. During cooling of the melt, the rates of nucleation 
and crystal growth in the material are greatly inhibited by the high viscosity 
of the melt. However if it is cooled rapidly to a temperature below the melting 
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point, the viscosity increases to such an extent that nucleation is totally 
prevented. The glassy state is preserved indefinitely. 

A fully crystalline sample could not be obtained despite using 
cooling rates as low as O.OlC/min, the addition of crystalline material to 
the melt to try to increase the rate of nucleation, and the maintenance of 
the molten material temperature very close to its melting point for 96 hrs. 
Measurements were finally carried out on non-crystalline void free samples 
fabricated into test slugs. 

As mentioned earlier ATJ graphite was used finally for preparation 
of the test sample in this case. After the molds had cooled the graphite 
pieces were broken away from the formed sample slug. The final machining 
and preparation of the samples had to be carried out with diamond tooling 
as for a normal glass. 

3.4 Experimental Procedure 

3.4.1 Differential Thermal Analysis 

A Netzsch Automatic Reocrding DTA instrument was used for these 
measurements . 

A small amount of the crystalline material was placed in a 
specially shaped nickel DDK measuring cup centered around a protected thermo- 
couple. A similar amount of powdered Kaolin as a reference material was 
similarly placed in a comparison measuring cup. The outputs of the two thermo- 
couples were connected differentially. A separate temperature measurement 
protected thermocouple was placed centrally between the two cups. The cups 
were placed at the center of a resistance heated temperature enclosure which 
could be heated or cooled at constant temperature rates. After allowing the 
system to come to equilibrium it was heated automatically at a constant rate 
of 5C/min and the output of the differential thermocouple recorded together 
with the temperature as measured by the measurement thermocouple until the 
measured temperature attained a value some ten to twenty degrees above the 
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the expected melting point or range. The whole system was cooled and a second 
sample measured using a temperature use of lOC/min. 

The melting point was determined from the intersection of lines 
drawn through the portions of the temperature time curve immediately below 
and above the peak. 

The system had been pre-calibrated with standard reference materials 

KC10., SiO„ and K„CrO, obtained from National Bureau of Standards. The measured 
4 2 2 4 

melting points were within + 0.5C of the reference melting points. 

3.4.2 Thermal Conductivity 

3. 4. 2.1 Solid Phase 


The basic method is shown schematically in Figure 3.1. For the 
measurements Pyroceram 9606, an inert ceramic glass, was chosen as the most 
suitable available reference material to cover all of the materials over the 
total temperature range. The thermal conductivity of Pyroceram 9606 is in 
the range 4 to 3 W/mK over the temperature range of interest and the salts 
were expected to be in a similar thermal conductivity range. 

In order to check the method and the apparatus a sample of Pyrex 
7740 glass was evaluated using the Pyroceram reference materials. This glass 
is another readily available reference material having a thermal conductivity 
in the range 1 to 1.5 W/mK. Results over the range 100 to 500C were within 
+ 1.5% of the accepted values for this material. These were considered most 
satisfactory as qualifying the methods and techniques to be employed. 

In all cases fine gauge chrome 1/alumel thermocouples in a high 
nickel alloy protective sheath having overall dimensions of 0.25 mm were fitted 
tightly into the fine grooves cut across the surfaces of the sample. The 
instrumented sample was sandwiched between two similar instrumented samples 
of the reference material and of similar dimensions. 
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STACK LOADING RODS 


INSULATION 



Comparative Method - Schematic Assembly 
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The composite sample was placed between two similar sized heater 
units and mounted on a fluid cooled heat sink under a uniform load applied 
at the top of the composite stack. A metal and ceramic heater guard tube 
some 75 mm diameter and 100 mm long which could be heated at various positions 
along its length and cooled at the bottom end was placed around the composite 
test stack such that the sample stack was positioned centrally with the guard 
tube. A further large metal fluid cooled shroud was placed around this stack 
and the interspaces around the sample and between the metal tubes filled with 
a thermal insulating powder which had been pre-dried under vacuum at elevated 
temperatures. A lid was fitted over the metal shroud and a large glass bell 
jar was placed around the whole assembly. The system was then evacuated and 
backfilled with pure dry argon. 

A steady temperature equilibrium was established in the system by 
means of adjustment of the power to the heaters in the sample stack and the 
rate of flow of cooling water through the heat sink. The temperature gradient 
along the length of the guard tube was matched approximately to that on the 
composite test stack by automatic control of the heater along the length of 
the guard tube. 

At equilibrium the temperatures in different sections of the stack 
were obtained from the various thermocouples in different sections of the stack. 
The thermal conductivity was derived from a knowledge of the heat flow as 
determined from the mean value calculated in the top and bottom reference 
materials, the temperature difference across the sample and the known dimen- 
sions as follows: 

(3.1) 


X = I & 

2 AT 


sample 


kAT 


Ax J 


top reference + 


kAT 


Ax J 


bottom reference 


where A = thermal conductivity 
Ay = thickness 
AT = temperature difference 

k = thermal conductivities of the reference at respective mean 
temperatures 


32 



Measurements were made at regular temperature intervals up to as 
close to the melting point as could be obtained without any part of the test 
sample being at a temperature above the measured melting point. Following 
these measurements a repeat determination was made at a lower temperature to 
check whether the sample had changed during the heating or for possible con- 
tamination of the thermocouples. 

A typical set of experimental data including temperature readings 
and calculation of results is given in Appendix B. 

3. 4. 2. 2 Molten Phase 

The method chosen was based on the comparative technique modified 
for the molten state. It is shown schematically in Figures 3.2 (a) and 3.2 (b) . 

Essentially the test sample is contained within a thin walled 
cavity in a larger piece of a reference material. The bottom solid section 
of the material acts as the lower portion of the comparative stack while an 
upper section of the same reference material fits tightly into the sample 
cavity to seal the test sample within the cavity. The upper reference material 
contained a vent hole to allow for expansion of the test sample on melting 
and subsequent increase in temperature. Sheathed protected fine gauge ther- 
mocouples were placed in the wall of the cavity, in the surfaces of the 
references touching the molten material, in holes along the length of the 
reference materials and in the sample itself. In all other respects the 
experimental techniques were the same as for the solid phase. 

In deriving the thermal conductivity of the molten sample the same 
principle was used as that for the solid material. However allowances had to 
be made for the heat flow through the walls of the sample container. This 
was found to be in the order of 20 to 25% of the total heat flow measured in 
terms of that in the upper and lower reference materials. 
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Molten Thermal Conductivity Test Stack Schematic Assembly 
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Numbers Indicate Positions of Thermocouples 



Figure 3.2 (b) 

Cross-Section of Pyroceram 9606 Test Stack 
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For lithium nitrate and one sample of the eutectic compound a test 
stack fabricated from Pyroceram 9606 was used. However the container cracked 
and no large enough piece of Pyroceram was available at the time from which 
a second test stack could be fabricated. In addition it had been found that 
lithium hydroxide, and the eutectic compound to a lesser extent, attacked 
pyroceram when they were molten and were assumed to be contaminated. 

A test cell of Inconel 600 was then fabricated for the measurements 
on the materials other than the lithium nitrate. The Inconel 600 was of much 
higher thermal conductivity (10 to 20 W/mK) than the Pyroceram 9606 but this 
factor was compensated by using a larger total cross-section and a thinner 
wall of the cavity. In this way the heat flow in the cavity wall was still 
kept to the order of 20% of the total heat flow. 

For the Pyroceram test stack a sample of calculated dimensions was 
cast and machined for insertion into the cavity. However for the Inconel test 
stack a sample could be cast directly in the cavity using the technique described 
earlier in Section 3.3. 

3.4.3 Heat Capacity and Latent Heat 

3. 4. 3.1 Adiabatic Calorimeter 


The apparatus is shown schematically in Figure 3.3. 

For the lithium nitrate and the eutectic material special sample 
containers were fabricated and calibrated empty separately. These were fabri- 
cated from 304 stainless steel in the form of cylinders approximately 50 mm 
long, 23 mm diameter with a 0.7 mm wall thickness with a screw cap. A resistance 
heater of sheathed wire was brazed all around the tube and on the top and 
bottom in a uniform configuration. The internal and external surfaces of the 
complete unit were gold plated and polished. Once a container had been cali- 
brated empty the specific heat of ice and water were measured as was the latent 
heat of fusion. Results well within + 1% of the accepted values for this 
material were obtained and the apparatus was then considered suitable for 
measurement on the salts. 
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LEGEND 

A SAMPLE CONTAINER - HEATER ASS'Y. 

B ADIABATIC JACKET 

C DIFFERENTIAL THERMOCOUPLE PROVIDINB 
JACKET CONTROL 8I0NAL 

0 THERMOCOUPLE FOR RECORDING SAMPLE TEMP. 


Figure 3.3 

Schematic Diagram of Dynatech QTA 
Quantitative Adiabatic Calorimeter 


37 




A sample was melted Into a container in the manner described 
earlier in order to maximize the amount of material being measured. The 
electrically heated gold-plated container was suspended inside a heavy nickel- 
plated and polished copper guard enclosure. The whole assembly was evacuated 
and then cooled at a uniform temperature below the lowest mean temperature at 
which data was required. When steady conditions were obtained, a controlled 
rate continuous power input was supplied to the heater on the sample container. 
By utilizing the output of a multi-j unction differential thermopile, the power 
to the guard heater was automatically controlled so that the temperature of 
the guard was equal to the temperature of the sample allowing negligible heat 
transfer from the sample to its surroundings. 

From observations of the power input to the heater and a continuous 
record of sample temperature variation with time, a record of the sample and 
container enthalpy change with temperature was obtained. 


The specific heat was calculated from: 


C 

P 


(mC ) 

E_ 


tot 


m 

s 


(mC ) 

R_ 


c 


( 3 . 2 ) 


where 


C 

P 


(mC ) 
p tot 


(mC ) 

P c 


m 

s 


= specific heat of the sample 

= instantaneous total enthalpy change with temper- 
ature from graphical record 
= instantaneous container enthalpy change with tem- 
perature from the calibration graphical record 
= mass of sample 


The latent heat of fusion was obtained from the direct measurements 
of the steady applied power to the sample and the time period at which the 
temperature remained constant. 


A typical set of data and its analysis is given in Appendix B. 
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3. 4. 3. 2 Drop Calorimetry 


The method used was based upon typical drop calorimetry procedure 
utilizing a calibrated copper sample receiver contained in a constant temper- 
ature water enclosure. 

Special sample containers were fabricated from 304 stainless steel 
in the form of thin wall cylinders some 75 mm long, 25 mm diameter and 0.7 mm 
wall thickness, with a threaded top cap . A container was calibrated at regular 
temperature intervals over the complete temperature range of the entire study 
and the other containers at individual points within the range. Slight dif- 
ferences in the overall enthalpy change for the different containers were 
adjusted from the knowledge of their masses. 

The specific heats of a sample of each of OFHC copper and of high 
purity alumina were measured with use of the containers as a final calibration 
of the system. 

For each test a sample was melted into a container as described 
earlier. The filled container was then attached to a 0.3 mm support wire 
of a nickel alloy and suspended vertically at the center of a three-zone 
controlled temperature furnace with the free end of the wire passing over a 
free moving pulley at the top of the furnace. Two thermocouples were suspended 
from the top of the furnace such that the junctions touched the sample near 
the top and bottom. The wire was so arranged that the sample was suspended 
vertically along the center line of the furnace and of such a length that 
when it came to rest at the end of its travel the sample rested upon the base 
of the receiver below it. Great care was taken to assure that the sample would 
fall quickly and freely without hitting the side of the furnace tube or the 
intermediate cylindrical protection tube and that it decelerated over the 
last 20 mm of its fall. 

The sample was allowed to attain the equilibirum temperature in 
a dry argon atmosphere for the order of 1 - 2 hours and regular readings of the 
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temperature measurement thermocouples taken. During the thirty minutes prior 
to a "drop" the temperature of the receiver was noted every 30 s. At equilibrium 
when the two thermocouples and the very small drift in temperature of the receiver 
had been noted, the sample was dropped quickly. As the sample dropped the 
radiation shields covering the intermediate zone moved sideways to allow the 
sample to fall and come to rest in the receiver. When the sample came to 
rest, these shields quickly returned to the original position in order to 
reduce any radiation transfer of heat from the furnace to the receiver or con- 
vective and radiant heat transfer from the receiver to the outside. 

The temperature of the copper receiver was taken regularly at 30 s 
intervals for the first thirty minutes following the drop followed by 60 s 
intervals for the next thirty to sixty minutes in order to determine the rise 
in temperature and the subsequent drift in temperature of the receiver. Fol- 
lowing a drop, the receiver system was allowed to come to equilibrium for the 
order of two hours. Two similar drops for each material at each temperature 
were performed. A typical data point is shown in Appendix B. 

Since in an adiabatic system AH r = -AH, (3.3) 

where AH^ = enthalpy change of the copper receiver due to sample drop at 25C 
-AH = total enthalpy change of the sample between the drop temperature 

and the final receiver temperature 
= mass of sample 

C = integrated average specific heat of sample over temperature range 
pa 

T to T, 
o f 

T q = temperature of sample in furnace 

T^ = temperature of sample and receiver at equilibrium temperature 

the calibration constant at each temperature for the copper receiver was calculated 
as : 


(mC ) 

P r 


AH 

r 



where (mC ) = receiver calibration 

P r 

T^ = initial receiver temperature 


(3.4) 
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The sample drops were made as indicated and the total enthalpy 
of the samples at each drop temperature was calculated as: 


AH = (mC ) (T - T ) (3.5) 

o p r i x - 

The total enthalpy of the samples at each drop temperature was plotted against 
the drop temperatures and the specific heat calculated as against the drop 
temperatures and the specific heat calculated as: 


mC 


dH 

ps [dt^ 


(3.6) 


where C = specific heat of sample at selected temperature 
= slope of enthalpy versus drop temperature plot 


m = mass of sample 


The heat of fusion was obtained in the following manner: 

1. Measuring the total enthalpy corrected to one receiver temper- 
ature, at two successive increasing temperatures one just 
below the melting point and one just above. The latent heat 
was derived in terms of the difference between the two values. 


2. Integrating the area under the curve of specific heat versus 
temperature over the temperature range immediately above and 
below the melting point. 


3.5 Linear Expansion and Density 
3.5.1 Solid Phase 


Originally the intention had been to use a combination of basic 
pycnometry and dilatometry for measuring density. However the pycnometer 
method utilizing the measurement of the change in volume of small chips of 
each material in an inert silicone calibration fluid could not be used due 
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to very slight reactions which took place when the salts were in contact with 
the fluid. In consequence the linear expansion of specimens of each type of 
salt, cut in mutually perpendicular directions from a large prepared solid 
sample, were measured in a dilatometer to enable the volume coefficient of 
each material to be determined. From the known mass at 20C of a given sample 
together with the known initial volume at 20C and its change with temperature 
the density at 20C and over the temperature range could be determined accurately. 


For the measurement of linear expansion the initial length of a 
sample was measured accurately with a micrometer. The sample together with a 
temperature measurement thermocouple was placed in a calibrated fused quartz 
push rod measuring system of a Netzsch Electronic Automatic Recording Dilato- 
meter. The system was placed at the center of a resistance heater environmental 
chamber and allowed to equilibrate. Power was then supplied to the heater in 
a regular manner such that the sample temperature increased at a constant slow 
rate of rise of 1 C/min up to as close to the melting point of the material 
as could be obtained before softening of the sample was noted. During the 
whole length of the experiment the continuous length and temperature changes 
of the sample were recorded. 


The coefficient of thermal expansion of the sample was obtained 

as follows: 


- AL 
LoAT 


(3.7) 


where a = the coefficient of linear expansion 
Lo = the initial length at 20C 

AL = the change in length for a particular temperature in- 
terval obtained from a curve of A1 versus temperature, 
allowance being made for the length change of the quartz 
system calibrated previously 
AT = the particular temperature interval 


The volume coefficient was obtained as the sum of the three respective linear 
coefficients. The system was calibrated using both a platinum tube and a fused 
quartz sample. 
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3,5.2 Molten Phase 


The basic buoyancy method was used to determine this property, 

A solid cylindrical piece of nickel some 6 mm diameter and 18 mm long weigh- 
ing approximately 40g at 20C was used as the suspended bob for the measurements. 
The expansion of a separate rod of the same nickel was measured in the Netzsch 
dilatometer in order to determine the change in volume of the metal bob at any 
particular temperature within the desired total range. 

A sample of the salt was melted into a large nickel container 
which could be maintained at any desired temperature level. The nickel bob 
was suspended from the arm of a balance into the molten salt. Shields and 
insulation were placed around the wire and between the furnace and balance to 
minimize radiation, convection and updraft effects. The mass of the bob was 
evaluated once the system had attained equilibrium. 

Measurements were made in this manner at successive increasing 
temperatures some 50 to 100C above the melting point. 

The density of the salt was determined as follows: 


P s = 


P Ni (m 20 


- m) 
T 


m 


20 


(3.8) 


density of salt 

density of nickel bob at test temperature (T) 
mass of nickel bob at 20C 
mass of nickel bob at test temperature (T) 

The system was calibrated using distilled water and a high temperature silicone 
oil. 


where p = 


m 


Ni 

20 

“T 


3 . 6 Results 
3.6.1 Melting Point 

The melting points of the particular samples of each salt tested 
are given in Table 3.3. 
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Table 3.3 


MELTING POINTS OF FOUR ENERGY STORAGE SALTS 
Melting Point, C 


Material 


Measured 

Published Handbook Values 

LiN0 3 (R) 
LiN0 3 (C) 



252 

252 

254 

63 LiOH - 37 

LiCl 

(R) 

262 


63 LiOH - 37 

LiCl 

(C) 

262 

262 

LiOH (R) 



470 

471 

LiOH (C) 



472 

Na 2 B 4 0 ? (R) 

Na 2 B 4 0 ? (C) 



737 

742 

741 


All of the above temperatures are very close to accepted values. This factor 
and the initial analyses confirm the identity of the materials. 

3.6.2 Thermal Conductivity 


The experimental results for the solid phase of the lithium salts 
are shown in Figure 3.4. Values of thermal conductivity obtained from smooth 
curves drawn through the experimental points are given in Table 3.4. These 
values are believed to be accurate to better than + 5%. 

The results for each indicate that the reagent and commercial 
grades are very similar in property when allowance is made for density and/or 
crystal form. It should be mentioned here that the reagent form of the LiN0 3 
was much less crystalline in appearance than the commercial grade whereas for 
the LiOH the reverse was the case. 

For the LiOH in particular the differences between the two samples 
are the largest. The difference in density is small and in the opposite direction 
to that which would be expected in order to explain the variation between the 
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Thermal Conductivity, X, W/mK 



Temperature, C 







Table 3.4 


THERMAL CONDUCTIVITY OF THREE ENERGY STORAGE SALTS IN THE SOLID PHASE 


Specimen 

Thermal Conductivity, W/mK at 

50 

100 

200 

225 

250 

300 

350 

400 

440C 

LiN0 3 (R) 

(p = 2139 kg/m 3 ) 

1.73 

1.53 

1.40 

1.38 






LiN0 3 (C) 

(p « 2215 kg/m 3 ) 

1.67 

1.46 

1.31 

1.29 






63 LiOH / 37 LiCl (R) 
(p = 1648 kg/m 3 ) 

1.33 

1.29 

1.21 

_ 

1.19 





63 LiOH / 37 LiCl (C) 
(p = 1635 kg/m 3 ) 

1.30 

1.27 

1.19 

_ 

1.17 





LiOH (R) 

(p = 1421 kg/m 3 ) 

2.50 

2.10 

1.65 

_ 

1.52 

1.40 

1.32 

1.27 

1.26 

LiOH (C) 

(p = 1404 kg/m ) 

3.06 

2.58 

2.00 

... . - . 

— 

1.89 

1.60 

1.49 

1.37 

1.35 















































two samples. However for solids in general the thermal conductivity is affected 
markedly by crystal size and, as mentioned earlier, the commercial grade sample 
appeared to consist of much larger crystals than the reagent sample. One 
other point, the commercial grade sample was prepared from the powder form 
whereas the reagent grade was from the crystalline form. This again could 
influence the thermal conductivity of a resultant test sample. 

The results for the borate material are quite different in nature 
to those for the lithium salt and are shown separately in Figure 3.5. Early 
attempts to prepare crystalline samples, coupled with results obtained in the 
specific heat measurement (to be discussed later) had indicated that it was 
not possible to fabricate crystalline solid samples and that anomalous enthalpy 
behavior occurred in the 500 to 600C regime. In consequence, measurements on 
the prepared samples in the solid state were carried out from approximately 
100C up to as close to the melting point as possible. 

It will be seen that above 500C there was an apparent large decrease 
in thermal conductivity followed by a further increase up to 700C. On dismantling 
it was found that a small amount of the sample in the form of beads had stuck to 
the sides of the lower test stack. This indicated that at some time during 
the measurements part of the sample had melted even though all of the sample 
was kept below 740C. As a consequence the thickness of the sample had changed 
but this was estimated to be 10 to 15%. An absolute measurement could not be 
made since the heat meter - sample - heat meter combination was a solidified 
entity and had to be broken apart and the surfaces of the heat meters remachined 
for the other measurements. 

The 15% maximum difference in thickness is not sufficient to increase 
the A value above 500C to fit to the line extrapolated through the experimental 
points below 500C. It would thus appear that the thermal conductivity of these 
samples does decrease significantly above 500C and that at some temperature 
between 500 and 700C limited melting of this particular form of sodium borate 
also takes place. 
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Results for the thermal conductivity of the molten phase for each 
material are summarized in Table 3.5. Some three data points for each material 
were obtained in one of the two cells described earlier. 

As for the solid phase there appears to be little or no difference 
between the reagent and commercial grades of each material. The results for 
the three lithium salts indicate an approximate halving of the thermal conductivity 
on melting. Such a decrease is in general agreement with results for other 
materials, including salts, when measured in the solid and molten phases. 

The discussion of the results for the sodium borate are necessarily 
obscured by the fact that changes in the thermal conductivity occur at temper- 
atures below the accepted melting point. These changes occur in fabricated 
samples which may be of a material which has supercooled to a large extent. 

However if one considers the results in the 500 to 600C range there is a decrease 
in the thermal conductivity of the measured sample of the order of 50% as for 
the other materials but it is believed that not all of the sample was molten 
at the temperature unless the very high viscosity of the material hold the 
sample in place in the test stack. Some material had escaped from the test 
stack during the measurements but it is not known at what temperature level 
this occurred. 

The thermal conductivity of the molten sodium borate measured in 
the molten cell at temperatures above 740C is similar to that measured in the 
solid samples at about 700C. They are believed to be representative of the 
molten phase. These values are only the order of 20% or so below that of the 
solid at 500C. However if one extrapolates the earlier result in the solid 
phase to 740C the reduction in value is some 30% or more. The whole question 
of the thermal conductivity of sodium borate remains obscure due to unknown 
factors which influence the form of the fabricated samples of interest par- 
ticularly below 700C and above 500C. 
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Table 3.5 


THERMAL CONDUCTIVITY OF FOUR ENERGY STORAGE SALTS 
IN THE MOLTEN PHASE 


Material 


Thermal Conductivity, W/mK 


260C 300C 


f- 5 ) 

\X / m.p. 


LiN0 3 (R) 

0.62 0.635 0.75 

0-45 

LiN0 3 (C) 

0.61 0.62 0.65 

0.49 


Material 

Thermal Conductivity, W/mK 

(“ ) 

\ X s /m.p. 

280C 330C 380C 

63LiOH/ 37L1C1 (R) 

0.70 0.735 0.78 

. 

0.57 

63L10H/ 37LiCl (C)‘ 

0.685 0.73 0.77 

0.57 


M n jil 

1 

Thermal Conductivity, W/mK 

iiaLci. laX 

490C 

550C 

600C 

LiOH (R) 

0.85 

0.87 

0.88 

LiOH (R) 

0.85 

0.88 

0.89 



Material 


Na 2 B 4°7 (R) 


Na 2 B 4° 7 (C) 


Thermal Conductivity, W/mK 


760C 800 C 


840 C 



50 


























3.6.3 Specific Heat and Latent Heat 


The results are summarized in Table 3.6. Enthalpy values for the 
LiNO^ and LiOH/LiCl eutectic are given in Figure 3.6 and for the other materials 
in Figure 3.7. 

The results for each material indicate that there were no differences 
between the reagent and commercial grades. The present values for the latent 
heat of the lithium salts are in very good agreement with previous measured or 
calculated values. The present specific heat results extend the temperature 
range for these materials over the complete range of use. Where previous 
values are available at lower temperatures the absolute values are in good 
agreement with the present ones. It is significant to note that the latent 
heat of fusion of LiOH is twice that of the other two lithium salts. 

Results for the Na^B^O^ samples do not show any indication of a 
significant latent heat at the nominal melting point. In fact, the results by 
drop calorimetry indicate a "heat of transformation" of the order of 1.5 x 10 
J/kg in the vicinity of 550C. It should be realized however that these 
measurements were carried out on pre-melted material and thus the samples were 
likely to be amorphous due to possible supercooling. The nature of the drop ’ 
calorimeter test, involving the rapid cooling of a sample from a fixed high 
temperature to a much lower one, would also give rise to further change of form 
or structure of a sample of this type of material. 

However examination of the DTA melting point curves carried out 
initially on the crystalline material had indicated that some reaction took 
place in the 540 - 560C region besides the nominal melting at 740C. It was 
decided to do a more careful DTA study on new samples of crystalline material 
over the lower temperature range and with the melt. 

A typical DTA curve is shown In Figure 3.8. It can be seen that 
a reaction does take place in the above lower temperature regime and that the 
material does melt at 740C. However, while these measurements are not accurately 
quantitative, it can also be seen that the heat of reaction involved at the 
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Tab le 3.6 


SPECIFIC HEAT, Cp, AND LATENT HEAT OF FUSION, h f 
OF FOUR ENERGY STORAGE SALTS 


C , J/kg K at h f 

P I 


Material 

160 

200 

240 

260 

280 

300 

340 

400 

440 

550 

600 

700 

750 

800C 

10~ 5 J/kg 

LiN0 3 (R) 

1540 

1640 

2000 

2050 

2170 

2270 









5.33 

LiN0 3 (C) 

1520 

1620 

1940 

2020 

2100 

2210 









5.28 

LiOH/LiCl (R) 

2100 

2240 

2370 

2440 

peak 

2650 

2800 








4.82 

LiOH/LiCl (C) 

2070 

2220 

2340 

2420 

peak 

2610 

2750 








4.88 

LiOH (R) 






2920 

- 

3150 

3380 

3940 

3660 




10.15 

LiOH (C) 






2730 

- 

3050 

3190 

4150 

3790 




11.10 

Na 2 B 4 0 7 (R) 









1300 

1400 

1720 

1755 

1775 

1795 

* 

Na 2 B 4°7 (C) 









1300 

1390 

1700 

1740 

1760 

1780 

* 


*Value could not be obtained at melting point for samples tested 



uiiLluxx* y , XU Lr, vr/ xvg 











TOTAL ENTHALPY, 10 h, J/kg 











*.551 / 
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Figure 3.8 DTA Curve for Sodium Diborate 
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lower temperature in some four times that at the melting point. If the approxi- 
mate value of 1.5 x 10 J /kg at 550C is assumed correct then the latent heat 
of fusion is thus only of the order of 0.4 x 10^ J /kg . 

After the original heating of the sample shown in Figure 3.8 had 
been completed it was allowed to cool very slowly to below 500C and then re- 
heated at the standard programmed rate. The resultant curve showed no indication 
of any significant changes from a straight line base curve in the 550C or the 
740C regimes either on cooling or reheating. 

Clearly the specific heat, latent heat, and the thermal conductivity 
of Na^B^O^ in the solid phase above 400C requires further study on well 
characterized samples if such could be made available. The present results 
however indicate that its use for energy storage purposes is very limited 
due to sample form limitations produced by different heating and cooling pro- 
cedures . 

3.6.4 Thermal Expansion and Density 

Summarized results of the thermal expansion of the solid materials 
are given in Tables 3.7 and 3.8 for the parallel and perpendicular to crystal 
growth directions respectively. As for the other properties there is little 
or no difference in expansion coefficient between the reagent and commercial 
grades. Similarly for the LiNO^ and the eutectic compound there is no significant 
difference in expansion behavior between the parallel and perpendicular directions. 
However for the LiOH and ^25^0^ there does appear to be a small but significant 
difference between the direction of each. In both cases the parallel direction 
is higher than the respective perpendicular direction. 

Since the differences were small, an average coefficient change in 
length was derived for each material and the volume coefficient was assumed to 
be three times this value. 
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Table 3.7 


THERMAL EXPANSION OF FOUR ENERGY STORAGE SALT MATERIALS ( 


Specimen 

4 

Thermal Expansion ^ (AL/Lo) X10 
Parallel to Crystal Growth 

20C 

50C 

100C 

150C 

20 0C 

225C 

300C 

400C 

450C 

500C 

600C 

lino 3 (R) 

- 

11.0 

30.0 

49.0 

68.0 

76.0 






lino 3 (C) 

- 

11.0 

31.0 

50.0 

70.0 

80.0 






63 L10H/37L1C1 (R) 

- 

10.0 

28.0 

46.0 

64.0 

73.0 






63L10H/ 37L1C1 (C) 

- 

8.0 

22.0 

37.0 

53.0 

62.0 





mi i 

LiOH (R) 

- 


21.0 


48.0 


74.0 

101.0 

114.0 


m 

LlOH (C) 

- 


20.0 


44.0 


69.0 

94.0 

106.0 


■m 

E 7a 

- 


6.0 


16.5 


29.4 

43.5 


57.4 

71.1 * 

»a 2 B 4 ° 7 

- 


6.0 


19.0 


34.0 

49.0 


64.0 

80.0 * 












IBI 


* Measurements not possible above approximately 600C due to softening of sample 



























Table 3.8 

THERMAL EXPANSION OF FOUR ENERGY STORAGE SALT MATERIALS 


Specimen 


LiN0 3 (R) 


LiN0 3 (C) 


63L10H/ 37L1C1 (R) 


63L10H/ 37L1C1 (C) 


LiOH (R) 


LiOH (C) 


Na 2 B 4 ° 7 (R) 


Na 2 B 4°7 (C) 


Thermal Expansion t (AL/Lo) X10 
Perpendicular to Crystal Growth 


34.0 


26.0 


13.6 


15.1 


c 

450C 

500C 

1 

600C 


56.0 

79.0 

48.0 

71.0 

26.5 

40.5 

28.0 

42.0 


54.9 

74.3* 

57.0 

77.0* 


Measurements not possible above approximately 600C due to softening of sample 





















































The density of each material at any temperature was calculated from 
the measured mean at 20C combined with the calculated volume at the particular 
temperature T as follows: 
m 

P T = -y— (3.9) 

T 

and V T = V 2Q (1 + 3a) (3.10) 

These calculated values of density are summarized in Table 3.9. 

Table 3.10 contains values of the measured densities of the salts 
in the molten condition at three temperatures above the melting point. Also 
included are values of the ratio of the density of the melt to density of the 
solid at the melting point. Since the sodium borate sample softened above 600C 
the value of the density of the solid at the melting point has been extrapolated 
for the results above 600C. 

Both LiNO^ and the LiOH/LiCl eutectic salts have a decrease in 
density of 23 to 28% but the change for LiOH is 2%. However there is an 
apparent increase in density of some 13% for the Na^^O^. This again is a 
further factor which makes the borate material unsuitable for energy storage 
purposes . 


In conclusion, after tests had been completed on the samples small 
pieces of each material were submitted for chemical analyses to see if any 
significant contamination had taken place. Some typical analyses are given in 
Appendix B. 


Little significant contamination was apparent. The LiNO^ and the 
eutectic materials contained some additional Si picked up from glass components 
used for part of" the preparation procedure or from the silica measuring system 
for expansion measurements. The other elements which were present in added 
amounts were iron and cobalt which were probably from the alloys used in the 
preparation procedure. However in all cases the total amounts were very much 
less than 0.1%. 
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Table 3.9 


DENSITY OF FOUR ENERGY STORAGE SALTS IN THE SOLID PHASE 


Specimen 

3 

Density, kg/m 

20C* 50C 

100C 

150C 

200C 

225C 

300C 

400C 

450C 

500C 

600C 

LiN0 3 (R) 

2180 2173 

2159 

2147 

2133 

2128 






LiN0 3 (C) 

2234 2227 

2215 

2204 

2191 

2184 






63Li0H/37LiCl (R) 

1589 1585 

1577 

1508 

1561 

1557 






63LiOH/37LiCl (C) 

1524 1520 

1515 

1508 

1501 

1498 






LiOH (R) 

1434 

1427 


1417 


1407 

1396 

1392 



LiOH (C) 

1376 

1370 


1362 


1352 

1343 

1330 



Na 2 B 4°7 (R) 

2354 

2351 


2343 


2334 

2325 


2315 

2304 

Na 2 B 4 °7 (C) 

2322 

2318 


2310 


2301 

2291 


2280 

2269 


* Density of actual piece of material tested, not 100% dense 




























Table 3.10 


DENSITY OF FOUR ENERGY STORAGE SALTS IN THE MOLTEN STATE 



Melting Point 

Material 

C 

LiN0 3 (R) 

252 

LiNO (C) 

252 

3 



3 

Density , kg/m 

zaoc 

310C 

34 OC 

1762 

1745 

1726 

1800 

1766 

1745 



Material 

Melting Point 
C 

63LiOH/73LiCl (R) 
63LiOH/37LiCl (C) 

262 

262 


Density , kg/m' 


290C 


1529 



350C 


1507 


/solid\ 

\Molten/ m -l 



Material 

Melting Point 


C 

LiOH (R) 

472 

LlOH (C) 

472 



500C 


1372 


Material 


Na 2 V 7 (R) 

Na 2 B 4 ° 7 (C) 


Melting Point 



Density, kg/m' 


800C 


2607 


Solid i 
Molten/ m * P ■ 




































Section 4 


ANALYSIS OF STORAGE SYSTEM 


4.1 Power Plant 

4.1.1 System Selection 

The impact of thermal energy storage was assessed for one of the four 
power plant systems. The high temperature gas reactor was not selected since 
the future of this approach is in doubt (Ref. 5). Any of the other three would 
be a satisfactory case for study. Thermal storage is likely to provide a 
greater economic benefit in conjunction with nuclear systems than for the 
fossil fuel fired system, since the heat production capital costs for the 
nuclear plants are a greater portion of the total cost. Also the projection 
shown in Figure 2.1 indicates that nuclear power will be the dominant part of 
at least one system in 1990. Therefore, the boiling water reactor was selected 
as the system for evaluation. 

4.1.2 Power Plant Specifications for Braun Sar 

Power plant data were obtained from the General Electric Company 
(Ref. 6). The cycle description is given on Figure 4.1. The system design 
calls for 1250 mw, being approximately 25% larger than the size specified in 
Table 2.1. This gives the design specification of Table 4.1. 

Table 4.1 

STORAGE SYSTEM DESIGN SPECIFICATION 



Absolute 

% of Daily Base Load 

Base Load Plant, MW 

1250 

100% 

Peaking Power Rate, MW g 

250 

20% 

Storage Capacity, MW e -hr 

480 

2% 


(Electric) 
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15.359.040* 

ll00P-420f-3»7 7h 


BRAUN SAR 
















4.1.3 Storage System Integration 


The cycle integration between the power plant and the storage 
system was investigated for two different modes of storage operation. In 
one case storage was assumed to occur over the period of five hours with a 
decrease in power plant output to supply the required heat input. In the 
other, there was a constant heat input to storage. 

Two cycles were investigated for the first case and are shown 

in Figures 4.2 and 4.3. The low pressure power turbine was operated at the 

design condition. Therefore, the power output of the low pressure turbine is 

the same, at 840 MW g . Thus the loss in equivalent power to supply the storage 

is 410 MW . In both cycles all bleeds and feed water heaters associated with 
e 

the low pressure turbine were maintained, and the steam to the storage system 
was condensed and returned to storage. The reactor output rate had to be in- 
creased since the enthalpy of the return stream was decreased. If the enthalpy 
of the return stream were maintained to keep reactor capacity level, the thermal 
energy to storage would equal the power loss from the high pressure turbine. 

In both cycles the steam flow to the low pressure turbine was 
throttled from reactor outlet. In Cycle 1 part of the remainder went to 
storage and part went through the system reheater to increase the temperatrue 
of the throttled steam to the design temperature. In Cycle 2, the entire 
remainder went to storage and the liquid condensate leaving the storage system 
was used to superheat the steam entering the low pressure turbine. Schematics 
for the two systems are shown on Figures 4.2 and 4.3. The performance is 
given in Table 4.2. 


64 



h = 1190.8 


Reactor 


= - j 

h = 534. 


• 
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Condenser 
& Feedwater 
Heating 



Low Pressure 
Base Load Turbine 


h = 534.5 


m = 


h ■ 


m = 10.52 


m 


= 3.2 


Figure 4.2 


weight flow 
lbs/hr x 10 

enthalpy 

Btu/lb 


Schematic and Heat Balance - Cycle 1 








Figure 4.3 

Schematic and Heat Balance - Cycle 2 



Table 4.2 


EFFECT OF THERMAL STORAGE ON CYCLE OPERATION 


Power Output, MW^ 

Reactor Heat Input, MW 

e 

Heat to Storage, MW 

e 

Equivalent Diurnal Electric 

Storage, 5 hrs Cycle, 25% 
Efficiency, MW^-hr 
Overall Efficiency 


Off Peak Storage 


Base Load 

Cycle 1 

Cycle 2 

1250 

840 

840 

3580 

3775 

4010 

0 

615 

890 


- 

770* 

1100* 

359 

.265 

.263 


* Cycle efficiency; // MW^—hr assumes about 65% of available thermal energy 
can produce useful output electrical energy (Ref. 1, Table 1). 

Both of these methods of integrating storage reduced the efficiency 
of the basic power system as well as imposing the conceptual inefficiency 
associated with storage. 

In the other mode of providing the storage the thermal output 
of the nuclear reactor was increased by 3% to supply the storage capacity. 

The flow rate corresponding to the second condition is given by: 

x 10 6 + 33,000 lb /hr (4.1) 


W 


pv 


3.412142 0 - 0.033 h + 16.2 
t gv 


h - h 
pv fw 


where W = vessel output steam flow (Ib/hr) 
pv 

= core thermal power (MW^) 

h = vessel output steam enthalpy (Btu/lb) 
pv 

h- = feedwater enthalpy (Btu/lb) 
fw 
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(Cycle data will be given in English units as supplied by General 
Electric. Only the final design data will be supplied in the metric system!) 

The schematic on Figure 4.4 shows the cycle operating conditions. 

Based on the discussion in Section 4.1.1 about the unlikelihood of 
reduction of power on nuclear systems, the cycles shown on Figures 4.2 and 4.3 
were eliminated. It appears more feasible to overrate the system by a small 
amount continually in order to provide a large peak load capacity. 

Thus the storage system will receive heat at the rate of approxi- 
mately 100 MW continually. It will supply heat at the rate of approximately 
800 MW during a peaking period of approximately two and a half hours. 

4 . 2 Evaluation of Materials 

The appropriate thermal and transport properties of the candidate 
materials for the four systems were determined and the results presented in 
Section 3. Three of the materials are viable candidates, as shown in Table 4.3. 

The fourth material, sodium borate, proved to be a completely 
unsatisfactory material for use in thermal storage, insofar as could be told 
by basic properties measurement of the pure material. When the material was 
initially heated from room temperature in the crystalline state there was a 
transition to an amorphous solid at approximately 820°K, with an associated 
change in enthalpy, and change in thermal conductivity. The heat of fusion 
associated with the actual heat atl010°Kwas much less than this. 

Once this initial transition occurred the material stayed in the 
amorphous state with no phase changes of any kind and thus with no heat of 
fusion. Some difficulty with fusion had been expected but not this much. 

Based upon the properties parameter LiOH and LiNO^ are approximately 
equal as candidates for thermal storage materials. Both are superior to the 
eutectic. Therefore LiNO^ is selected for system design at the specified 
temperature level. 
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\D 


m 

h 



flow rate 
lbs/hr x 10 

enthalpy 

Btu/lb 


Figure 4.4 


Constant Heat Supply System 






Table A. 3 


SUMMARY OF MATERIAL PROPERTIES 


Material 

LiN0 3 

LiOH/LiCl 

LiOH 

Melt Temperature °K 

527 

535 

7A3 

Density 

Test 

2.310 

1.900 

1.A15 

P,kg/1 

Postulated 

(2. A) 

(1.6A) 

(1.A6) 

Heat of Fusion 

Test 

530 

A85 

1080 

h p , J/kg 

Postulated 

(379) 

(437) 

(930) 

Thermal Conductivity 
X, W/mK 

1.37 

1.10 

1.20 

Volume Product 

Test 

122A 

922 

1A86 

ph p , J/l 

Postulated 

(910) 

(717) 

(1358) 

Properties Parameter 
Ahpp 

A1 

32 

A2 
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4.3 Storage System Design 

4.3.1 Specification 

The storage system design will be based on Figure 4.4. The storage 

system will continually receive heat from the power system in the form of a 

steam bleed. During the peaking period of 2.7 hours the system will receive 
water from the peaking power system. System design operating conditions are 
shown on Figure 4.5. 

Steam will flow through a recuperator containing static heat of 
fusion material. During heating the steam side is maintained at reactor pressure 
and the heat of fusion material is liquefied. During the power cycle the steam 
pressure is dropped to the low pressure turbine pressure and the heat of fusion 

material is solidified, while supplying heat to the low pressure system. The 

bypass steam flow is maintained at all times, being throttled to mix with the 
peaking steam flow during peaking operation. Make up condensate from the 
peaking cycle must be pumped into the reactor circuit to maintain liquid quantity. 

With this design configuration the specification for the heat 
transfer system is given in Table 4.4. Table 4.4 shows the peaking period is 
controlling and decides the heat exchanger size. It also shows that the design 
is well balanced, unless there is natural convection in the liquid forming the 
liquid side heat transfer. 

4.3.2 System Design 

The heat exchanger design can be based on material inside the tubes, 
in a flat plate heat exchanger, or outside of the tubes. 

4. 3. 2.1 Storage Tubes Design 

With an idealized design the tubes will be sized to contain the precise 
amount of volume when liquefied. This specified the ideal radius ratio at the end 
of the solidification. Using this ratio the diameter and length of tubes could be 
determined. 
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m = 0.55 x 10^ lbm/hr 


Reactor 

Steam 

Flow 


Peaking 

Power 

Flow 


h = 1190.8 Btu/lb 

T = 540. 3F 
s 

T = 540. 3F 

— VV^A/^/v^AAAAAA/V^AA/^A/VV^AAAAA/ — 

s 

p = 965 psia j q = 104 MW 


Storage Material 


T = 527K = 489F 
tn 


Storage = 2215 MW-hr 


q = 833 MW 


m = 3.275 x 10 lbm/hr i 

h = 275.5 Btu/lb 

i 



— AAAAA^VvAAaAAAA/WVVVVVXAAAaA/-^— 


T = 370. 3F 

T = 310F 

S 


p = 174.2 psia 


Period of flow = 2.7 hours 


Figure 4.5 

Storage System Design Specifications 


= 540.1 Btu/lb 

■ T = 508F 

s 


h = 1252.1 Btu/lb 
T = 464F 


Table 4.4 


STORAGE DESIGN CONDITIONS 


Capacity 2215 MW-hr 

Heat of fusion, J/kg 0.53 x 10^ 

Period 

Storage 

Peaking 

Heat load, MW 

104 

853 

Temperature difference, °K 

28 

66 

Time, hr 

21.3 

2.7 

3 

Density, kg/m 

1800 

2310 

Thermal Conductivity, W/m °K 

0.62 

1.37 

Area, m^ 

.224 
x 10 6 

.243 
x 10 6 

Design melt thickness, m 

.0374 

.0268 
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r 

o 



1.8 

2.31 


= 0.47 


(4.2) 



From Equation A. 10: 


r 

o 



2AATt 


<r i /r o ) 


7 

| i n 

[r /r 



l ° lj 


V 2 x 1.37 x (530 - 461) 2.7 x 3600 

.53 x 10 6 x 2400(1 - . 47 2 ) ln(l/ .47) 

= .05 m ( ~ 2 inch) 


Lithium nitrate did not demonstrate any of the corrosion problems 
associated with lithium hydroxide and lithium chloride. Therefore, stainless 
steel tubing rather than inconel is specified with an eighth inch wall (0.32 cm) 
leaving an inside diameter of 1.75 in, with a ten percent margin of error on 
the design diameter to allow for shell side thermal resistance. 


The total amount of tubing required is specified by the liquid 
volume. The total volume required is given by: 

v _ Q _ 833 x 10 6 x 2.7 x . 3 600 ^ (4.3) 

h F P £ .53 x 10 6 x 1800 

- 9550 m 3 


The total length of heat of fusion tubing is: 


V 

A 


9550 
TT ( . 05) 2 


= 1.22 x 10 6 


m 


(4.4) 
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The shell volume for center to center triangular spacing of 1.125 times the 
outer diameter is : 


V 

s 



(4.5) 


V = 9550 x 
s 


1.125 


.875 


6 




= 17407 m 3 

With a tube length of 50 meters, the pressure vessel diameter becomes: 


/ 17407 X 4 
" K 50 tt 

= 21 meters 


(4.6) 


This is not feasible even with a number of vessels. 

4.3.2. 2 Design with Steam in Tubes 

With steam inside the tubes there is no limitation on the radius 

ratio in Equation A. 10. As smaller tubes are selected the area requirement 

for heat transfer becomes smaller and thus less expensive. Thus an arbitrary 

design radius selection of .0125 m (1 inch O.D.) was made with a .0025m wall. 

For design purposes the steam side heat transfer coefficient was assumed to 
2 2 

be 2500 W/m °K (~500 Btu/hr ft F) with a stainless steel thermal conductivity 
of 19 W/m°K. This reduced the effective driving temperature difference for 
solidification in Equation A. 10 by the temperature difference required to 
transfer 833 MW through these two resistances at the given area. Equation 
A. 10 was solved by trial and error for the radius ratio of the solidifying 
interface at the end of the peaking period. Then the area was determined from 
Equation A. 9 and used to define the difference more precisely. The final 
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resulting design was: 


Tube outer diameter 

Radius ratio at maximum 
solid layer diameter 

Tube area 

Tubing length 


0.025 m 
3 

133,600 m 2 
1.7 x 10 6 m 


The recommended tubing area would be triangular staggered with the minimum 
tube to tube spacing being determined by liquid volume corresponding to the 
solid mass frozen. If at least this much space is not provided there would 
be a possibility of liquid being trapped within a solid mass and imposing an 
excess stress on the tubes. Thus the minimum tube to tube spacing (D^/D^) 
is specified by: 



= 3.21 (use 3.25) 


(4.7) 


With this spacing the total volume required for the tube bank is: 

V-r s 2 ^L (4.8) 

= (3.25 x .0125) 2 1.7 x 10 6 

= 9720 m 3 
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In addition, an excess liquid level must be provided to cover the top of the 
tubes at all times. This is equal to the volume difference between the liquid 
and solid volumes of the active storage volume: 



s tt(.0125) (3 “ 1) 1.7 x 10 1300 *"* ^ 

- 1890 m 2 


The total mass of material required is: 

" - <V t ■ * r o 2uc i 

o C. 

= 1890 + 9720 - tt(. 0125) 1.7 x 10 x 1800 

, - 

= 19.2 x 10 6 kg 
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Appendix A 

DERIVATION OF AREA CRITERION 


The solidifcation rate of a pure substance on a solid interface, 
when the interface is maintained at a constant temperature, depends upon the 
rate of heat transfer through the solid layer. The liquid can be treated as 
a well stirred fluid, with relatively little heat transfer to the superheated 
liquid. The sensible heat transfered to the subcooled solid is small as com- 
pared to the heat of fusion. As a first approximation in the following analysis 
it is neglected. A more precise, slightly conservative, approach would be to 
consider the sensible heat in the solid layer as being generated at the solid- 
ification face. This quantity, which closely approximates one-half the product 
of the solid specific heat and the temperature difference between melt and 
surface, can be added to the effective heat of fusion. 

The heat flux at the solidifying face is: 

q/A = ph p (A. 1) 


With the assumptions specified above: 


q/A = ^ (T - T ) 
x m s 


(A. 2) 


Combining A.l and A. 2 gives: 


. (T - T ) 
dx a m s 


dt x h p p 


(A. 3) 


Integrations over a period of time (t) give a thickness: 


A(T - T ) 


xdx = 


m 


V 


dt 


(A. 4) 


A(T - T ) 
m s 

h F P 


(t l ' t o ) 


(A. 5) 


AO 



Thus the solid thickness as a function of time for these con- 


ditions is: 


A 


J A(T m- T s )(t l- ^ 

M v 


(A.6) 


The total amount of heat transfered, as determined from the total 
thickness solidified is: 


Q 

A 


= h pP A 


= /Xk 


F p(T m ~ 


V (t i - 


t 0 ) 


(A. 7) 


Thus the total amount of heat transfered through a solidifying 
layer is proportional to the product of properties grouping, /Xh p p , and 
system design conditions. Obviously this approach to system analysis is 
simplistic since there will be an interaction between the rate of heat transfer 
at any instant in time and the interface surface temperature. However, it 
does separate the effect of material properties for the purpose of a first cut 
at material selection. 

In the analysis of solidification on or in a tube, an explicit 
solution cannot be developed. If the limiting case is at the completion of 
solidification, the heat transfer rate is: 


5. = 


AAT 


A r ln(r /r) 
o o 


(A. 8) 


The total heat transfer is: 

2 - V (r Q 2 - f2) 
A 2r 


(A. 9) 


This produces an implicit function for the radius ratio: 


1 - (r/r )‘ 
o 


- y / \ i 2 A ATt 

ln<r o /r) 2“ 

r o V 


(A. 10) 


A1 



The area is inversely proportional to the function of the thermal 
conductivity: 


qr ln(r / r ) 
o o 1 

XAT 


(A. 11) 


Thus, although the properties enter in the same manner as for the rectilinear 
use, it can not be demonstrated that the simple form of the properties para- 
meter applies. 
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APPENDIX B 


Bl-2 

B3-4 

B5-6 

B7-10 


Test Data and Analysis 


Typical Data Point For Thermal Conductivity 
Typical Adiabatic Calorimeter Analysis 
Typical Drop Calorimeter Analysis 
Typical Chemical Analysis of Materials After 
Measurement 
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DYNATECH CORPORATION 


DATA REDUCTION SHEET 


Thermal Conductivity by Comparative Method 

Temperature Point: g<=}o C- Project: A/AS-S2-' 

Sample Designation: ^ T<L . 3 A Thickness, X Q : ( )-( 

Bottom Heat Meter Material: Thickness, X BHM : ( ) - ( )= (//ff,)cm 

Top Heat Meter Material: Thickness, X THM : ( ) - ( ) = (//^)cm 


Thermocouple Material: 


Avg. of (how many) readings: 3 Data of (date): 



Temperature Drops k Heat Meters, watts / m deg C 

BHM: BHM: £.Zf @3* Z. °C 


Sample: 44 tOC> 

THM: H THM: £ Z Z @ «^>/ °C 



Energy Spread 


Un 


/. 





















DYNATECH CORPORATION THERMAL CONDUCTIVITY DATA SHEET 



Additional Information: 


Page 


of 


mm 















DYNATECH CORPORATION 


QTA DATA REDUCTION 


//4/7& 


Sample: /./A/0 


Power: .535 


V 


Data of: z/zo/70 Project: 5% ~ / 


r> i" _ i ' Weight before:, M?7f kg: after: .60/83 k 8: (.OC&Z\ 


watts x 60 sec/mln = 2 , / watt - sec/mln 



-145.4 


-124.4 


-105.4 


- 87.9 


- 71.4 


- 55.9 


- 41.1 


- 26.9 


- 13.2 


0 


12.8 


25.3 


37.4 


49.2 


60.8 


72.1 


83.2 


94.0 


104.7 


115.3 


125.6 









a 

.04855 

1 


MfWfi 

c 

.0493 

Kim 


39. V 

0 

.0499 

wnm 

'! 

Wm 

r 

.0506 

mmm 

mm 

1) 

35.31 

0.0512 

II 

35. 67 


)3.17 


13.79 | M3 



KflHBBSUIEffTMKlyjJfc 


185.1 




8.50 


0.0523 












































































































































DYNATECH CORPORATION 


Date: Z //?/?£ Data of: 2-//3/Z ' 


QTA DATA REDUCTION 


Power: 


watts x 60 



dV 

dV 

dt 

dT 

mv 

mv 

min 

degC 

.6%6 

0.0529 


dee C 


min 


watt s 


min 



degC 










deg C I deg C [ kg deg C 


!C30 




V. 55 I 


H5 































































































SAMPLE # a. t, 0 » -_ a%7 , t 


■ r~ 


3 




DROP CALORIMETER 

7 >* yju± 


PROJECT: /i/ficr 
DATE: 


IE 


DIFF. 


OIL 















10 X 10 TO THE CENTIMETER IB X 2b CM 
KEUFFEL & ESSER CO in u b a 




/ 


_nOn 


JARRELL-ASH DIVISION 

ObO Fisher Scientific Company 


OX^rltfiratr of Attalgata 


TQ; Dynatech R/D Company 


DATE RECEIVED- 


99 Erie Street 


Cambridge, MA 02139 


DATE REPORTED:. 
ORDER NO.: 


Attn: Mr. A. Desjarlais 


SAMPLE DESCRIPTION:. 
INSTRUMENTATION: 


Li NO 3 


3.4 Meter Ebert Mark IV Spectrograph 


8/2/76 


8/6/76 


13597 



TL 21*103 


590 LINCOLN ST., WALTHAM, MASS. 02154 


litho in U.S> 

B7 















































































































JARRELL-ASH DIVISION 

Fisher Scientific Company 


(Eerttftrate of Analysts 


SAMPLE DESCRIPTION: 
INSTRUMENTATION. _ 


Dynatech R/D Company DATE RECEIVED . 

99 Erie Street DATE REPORTED: 

Cambridge, MA 02139 ORDER NO.: 

Attn: Mr. A. Desjarlais 

SCRIPTION: Eutectic 

itation 3.4 Meter Ebert Mark IV Spectrograph 


8/2/76 

8/6/76 

13597 



500 


1 


100 


200 500 


5 

6 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 



ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 


5 

6 . 

ND 

ND 

ND 

ND 

ND 

ND 

5 

5 



REMARKS: 


500 


10 


5 = Reagent; 6 = Commercial; Results in ppm 
‘except where percent is indicated. 


Lu 


Hf | ND 


ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 



KEY: 



ND - Not Detected 

T 

.01 -.1% 

VVFT < .0001% 

L 

. 1 - 1% 

VFT .0001%- .001% 

M 

1% -10% 

FT .001% - .01% 

H 

>10% 


STAFF ANALYST 


SUPERVISOR. TESTING LABS 


1 ! 21-103 


590 LINCOLN ST., WALTHAM, MASS. 02154 


litfto m U.S-^- 

B8 




















































































































£ JARRELL-ASH DIVISION 

JcP Fisher Scientific Company 


(Efrtifiratc of Analysts 


DATE RECEIVED 
DATE REPORTED:. 
ORDER NO.: 


SAMPLE DESCRIPTION: Li OH 


INSTRUMENTATION: 


1.4 Meter Ebert Mark IV Spectrograph 




= Reagent 


ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 



ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 


= Commercial 



Results in ppm. 


ND 

ND 

ND | 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 





KEY: 

NO - Not Delected T .01 -.1% 
WFf<.000l% L .1-1% 

VFT .0001%- .001% M l%-IO% 

FT .001% - .01% H >10% 


STAFF ANALYST 


SUPERVISOR. TESTING LABS 


TL 21-103 


590 LINCOLN ST., WALTHAM, MASS. 02154 


Litho in U.S>- 
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JARRELL-ASH DIVISION 

Fisher Scientific Company 


(tettftrate of Attalgais 


Dvnatech R/D Company 


DATE RECEIVED. _ 
DATE REPORTED.. 
ORDER NO.: 


SAMPLE DESCRIPTION: 
INSTRUMENTATION: 


Cambridge. MA 02139 ORDER NO 

Attn; Mr. A. Desiarlais 

SCRIPTION : Na 2 B 4 °7 Reagent 

ITATION ■ 3.4 Meter Ebert Mark IV Spectrograph 


8/2/76 

8/6/76 

13597 



REMARKS: 


Results in ppm except when percent is indicated 


VVFT < .0001% L .1- 1% J 

VFT .0001%- .001% M 1% -10% j 

FT .001% - .01% H >10% ' 


STAFF ANALYST 


SUPERVISOR. TESTING LABS 


’l ?}-103 


590 LINCOLN ST., WALTHAM, MASS. 02154 


Litho in O S A 

BIO 


























































